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ABSTRACT 


The combustion of a porous gravhite fiber olate is analyzed. 
Pearanslent One-dimensional, mathematical model is utilized to 
conduct numerous computer tests. The model simulates the ther- 
mal response of a porous graphite fiber plate and it is based 
On apnropriate energy equations. Two of the energy equations 
are based on energy balances conducted on both the porous struc- 
ture and on the air within the graphite fibers. The third 
energy equation is based on a mass balance performed on a 
differential volume of porous medium. 

A test plan is formulated and computer test runs are con- 
M@metea tO investigate the effects of plate thickness, filament 
diameter, permeability, porosity and exterior wind velocity on 
Brew thermal response of the porous graphite fiber structure. 


Mie test results are analyzed and conclusions are presented. 
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i eercesuGrroN TO" THE PROBLEM 


The advent of advanced fiber-reinforced composite materials 
mas Spurred substantial research efforts across a broad spec- 
trum. Whole organizations have been formed to analyze, design, 
and fabricate components made of composite materials. 

mEpGous Composites consist of long fibers contained within 
@ Matrix material. Long fibers are inherently much stiffer 
ma stronger than the same material in bulk form. Fiber pro- 
Berties differ from those of the bulk form due to the more 
perfect structure of the fiber. Additionally, fewer internal 
defects are present in fibers than in bulk material. Strengths 
and stiffnesses of a few selected fiber material are presented 


in Table I below: 


Wrobel: FLBER PROPERTIES 


TENSILE TENSILE 
DENSITY,o STRENGTH,S S/o STIFFNESS, E E/p 
FIBER Maen | wie? lb/in-| [10°in] [10*1b/in*] [10’in] 
Aluminum .097 90 9 LO sie iit 
Titanium 7 0 280 16 hee 10 
Steel me 2 600 oil Se) ial 
Boron .093 500 54 60.0 65 
Graphite me Sil 250) 49 S07 -. 0 a2 


Boron and graphite fibers are frequently selected for com- 


DPO-wtes that will be uthlized as aircraft structural material. 





The matrix material in such applications is frequently a plas- 
tic such as epoxy or polyimide. These advanced composites are 
characterized by properties that are superior to those of ma- 
terials that have been previously utilized for aircraft struc- 
tures (1.e. aluminum and titanium). Two major advantages are 
mmproved strength and stiffness, when compared with other 
materials on a unit weight basis. 

Advanced composite materials have been designed that are 
three times as strong as aluminum, yet weigh only sixtv percent 
as much [Ref. 1]. It is this favorable strength to weight 
ratio that makes composite materials particularly attractive 
to the aerospace industry. 

Currentlv, almost every aerospace company is developing 
Meoducts made with fiber-reinforced composite materials. 
Replacement components are frequently made with composite 
Materials. Examples of this application are the boron/epoxy 
fuselage section and horizontal tail on the Seneral Dynamics 
F-lll. Several companies are also introducing production 
components made from composite materials on their new aircraft. 
An example of this application is the F/A 18 Hornet Strike 
Fighter currently being built by McDonnell Douglas. The F/A 18 
airframe is a balance of conventional materials and graphite 
epoxy composite. The wing skins, trailing edge flaps, 
Stabilators, vertical tails and rudders, speed brake, and many 
access doors are made from composites. This balanced design 
has resulted in appreciable weight savings and increased 


aireraet performance. 








Weight Savings are of critical importance in military air- 
craft applications. Weight reductions result in the following 
desirable features: 

(1) Increased payload 

(2) Improved operating efficiency 

(3) Extended range 

[t is these attractive features that have created a great 
meal of interest within the U. S. Navy. Lighter weight air- 
craft with increased payload are ideal for carrier operations. 
Additionally, composites materials have improved the feasibility 
of VSTOL aircraft, which could lead to smaller, less expensive 
mimecraft carriers. 

ime Imeereduetion of aircraft structural components made 
from composite materials and the potential of aircraft made 
entirely from composites has spvurred a great deal of Navy 
Sponsored research activity. Research is being conducted to 
gain insight into the following features of composite materials: 

(1) Thermal properties 

fe) Strenctheand stiffness 

(3) Fatigue life 

@) Gombustion characteristics 

The combustion characteristics of composite materials and 
their thermal response in a fire environment are of particular 
Semeermm to the Navy. The potential for aircraft carrier flight 
deck pool fires is quite high, due to the presence of various 
volatile fuels and ordnance on the flight deck during flight 


operations. 
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Pees nesponse Of metallic materials to aircraft carrier 
flight deck pool fires and the reauired fire-fighting tech- 
meaues are well-known, but little data is available concerning 
tne combustion characteristics of composite materials. In 
meaer tO gain an insight into these characteristics an initial 
series of tests was conducted at the Naval Weapons Center, 

China Lake [Ref. 2]. The objective of these tests was to de- 
termine the response of graphite-epoxy composite materials to 
carrier flight deck pool fires and compare them with equivalent 
Strength aluminum. 

For these tests, a composite test specimen was mounted in 
such a manner as to simulate the heating condition of a complete 
Wing subjected to a large pool fire. The various test specimens 
were then subjected to a steady state thermal environment. | 
Specimen temperature data were recorded by thermocouple instru- 
mentation until all internal reactions had ceased. 

The combustion process was observed to proceed in two 
Sistinct phases: 

mm) Epoxy combustion 

(2) Graphite combustion 
mee tie completion of the first stage, the epoxy was totally 
consumed, exposing the porous graphite fiber structure (epoxy 
ignites and hurns at a much lower temperature than graphite). 

In some cases, the remaining graphite fibers exhibited 
mimeracteristics of self-sustained combustion and a tendency 
to smolder, while in others the fibers would cool, indicating 


@m inherent resistance to fire burn-through. 


Jal 








These initial tests indicated that the thermal response of 
mmemmoosite plate 1s not constant, but is a function of the 
MariOus parameters that characterize the composite material 
(i.e. filament diameter, ply thickness, plate thickness), the 
exterior wind velocity, and the thermal environment. The 
number of parameters involved in the problem indicated that 
an exhaustive experimental program designed to further inves- 
tigate the thermal response of a composite plate would not be 
economically feasible. Therefore, a comnuter model appeared 


to be a desirable means for conducting numerous tests. 


a2 





mie Shoe WON Or MATHEMATICAL MODEL 


In order to conduct a thorough computer analysis of the 
combustion characteristics of a graphite-epoxy composite plate, 
an appropriate mathematical model was required. A suitable 
one-dimensional model was previously formulated by Vatikiotis 
(Ref. 3]. The model was to simulate the thermal response, in 
@m@e-dimension, of the graphite fiber mat after epoxy consumption. 
The following discussion summarizes the assumptions, governing 
equations, boundary conditions,and method of numerical solu- 


mon utilized in the model. 


Me AVERAGE PORE VELOCITY 

In formulating the model it was assumed that an exterior 
wind velocity would always be present, thereby providing air 
flow over one surface of the composite plate. This assumption 
1s quite realistic when one considers the relative wind velo- 
memes consistently present on aircraft carrier flight decks. 
The air flow produces a pressure differential across the com- 
posite plate. This pressure differential induces a convective 
air flow through the porous graphite structure governed by 


Darcy's Law, 


u, -2¢ COE) 


: ; ; . ’ dP 
where u_ is pore velocity, uw is the dynamic viscosity, / ax 
[etme pressure gradient, and m is the specific permeability 


of the medium. 


AES: 





As a result of extensive experimental work, it has been 
suggested that farcy's Law is valid only in those instances 
where Re; < 1 [Ref. 4]. Therefore, internal Reynolds number 
Was Calculated in all test cases conducted. In all instances 
Re; mes On the order of 1 or less. 

Bene vVelocitv is a critical factor in the development of 
this mathematical model. This air flow enhances internal con- 
vection heat transfer, providing a means for extinction and 
eventual cooling of the graphite fibers. The air flow is also 
a source of oxygen for self-sustaining combustion, providing 
a means for heat generation within the porous structure. If 


the heat transfer, at a point, due to the air flow is greater 


meen the heat generation, then cooling will occur at the point. 


If heat generation is greater than heat transfer at a point, 
then combustion will proceed at this location. 

Equation (II.1) indicates that pore velocity OF? nis el 
function of permeability (m). As a result of the critical 
relationship between pore velocity and thermal response, it 
was immediately assumed that permeability would be a critical 
factor in the formulation of the mathematical model. 

Specific permeability is indicative of the ability of a 
porous medium to allow flow. It is dependent on geometry and 


on the properties of the medium itself. Several empirical 


models have been developed to obtained values for permeability. 


Vatikiotis selected a serial tvpe representation proposed by 


Bemeidegger {Ref. 5] in formulating this mathematical model: 
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m= 3g p(S/t)? (II. 2) 
where the porosity, p, is defined as, 
7 2 2 


mes) the filament diameter and D is the ply thickness (see 
Meeure 1). The average pore diameter, 6, was calculated using 


am @malogy to the hydraulic radius commonly utilized in fluids 


dyanamics, 
eeeenicd 
SP ARCTEN(S) tte, 


The tortuosity, T , is a kinematical property equal to the 
average length of the flow path of a fluid particle from one 
Side of a porous medium to the other. It is also a geometrical 
property dependent on the ratio d/D. 

The above discussion indicates that the pore velocity 


utilized in the mathematical model is effected by the following 


parameters: 
(1) Plate thickness (D) 
(2) Filament diameter (d) 
(3) Permeability (m) 
wy Porosity (Dd) 
Poeeteester ior velocity (a 


B. GOVERNING EQUATIONS 
After an expression for average pore velocity had been 
derived, appropriate energy equations were developed. 


Vatikiotis [Ref. 3] performed energy balances on both the 


i 








Gemous Solid and on the air within the graphite fiber structure 
in order to develop the proper equations. The energy balances 
Membaged the control volume approach. This approach is sum- 
Mmericed as follows: 


Heat into + Heat Generation = Heat out + Increase in Internal 
dv within dv of dv Bile Gey 


where dv = dxdvdz. 

By utilizing the porous medium as a control volume, the 
fiber heat transfer equation was formulated. Conduction, 
convection, radiation, change in internal energy, and heat 
generation rate per unit volume were all considered in the 
formulation of this equation. 

In developing the fiber heat transfer equation, Vatikiotis 
[Ref. 3] Limited the model to the kinetic regime of combustion 
and assumed that the dominant chemical reaction was 

co tate > C5 
This idealized treatment of the actual reaction is a good 
approximation during the initial combustion of the graphite 
fibers. 

An internal flow heat transfer equation was obtained by 
performing an energy balance on a differential volume of air. 
Heat conduction through the air, convection from the fibers 
Gemtcme dir, and energy transport due to air flow were all 
considered in the formulation of this equation. 

In developing the internal flow heat transfer equation, 
Vatikiotis [Ref. 3] utilized the perfect gas assumption for 


air in order to calculate enthalpy as follows: 
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dh = CaT. fo lees 50) 


where oe Pomeeopecitic meat of air at constant pressure. 

The effects of the combustion by-products on air were 
assumed to be negligible allowing the mass flow to be repre- 
Sented as, 


= 2g be ght. ©) 


Equations II.5 and II.6 were then utilized to determine 
mime energy transport due to the air flow, mh. 

In both the fiber heat transfer equation and the internal 
flow heat transfer equation, thermal conduction was determined 


by Fourier's Law of heat conduction in the following form: 


eT 
Geond = pas dvdz (Ei) 


Similarly, convection heat transfer was determined in both 


equations by application of the following expression: 


q nee) i ti dA (ites) 


conv im < 
where T_ is the temperature of the graphite, T,_ is the tempera- 
g i‘ 2 

ture of the air, dA is the surface area of the graphite in the 
differential volume, and h; fine imeernal camvection heat 
Smanster coefficient. An empirical expression developed by 
Yoshida, Ramaswami, and Hougen [Ref. 6] was used to determine 
Dus. 

1 

The final field equation formed was the oxygen transport 

equation. In order to formulate this equation a mass balance 


was performed on a differential volume of porous medium. The 


mass balance was formulated in the following manner: 


ily 





_ 


ae 


Oxygen into = Oxygen out + Oxygen Consumption + Oxygen 
dv of dv Accumulation 


Two mass transport mechanisms were considered in formulat- 
ing this equation: 

fei) Mp = molecumar diffusion 

fig, ) Me MecnwVc@mEl~c transport due to mass flow 


Molecular diffusion was determined by Fick's Law, 
m, = Be2 dydz ere) 


where B is the diffusivity of oxygen into the porous medium. 
The convective transport term was determined by 


m = pee ce Cert t)) 


Oxygen consumption rate was determined by multiplying 
@eeeeton Tate, r({T_, $), by the stoichiometric ratio (1/f) of 
o 


the reaction. The reaction rate is of the Arrhenius type, 
feel yd) exp(-E/RT) 


where ~ a Pomeneomeiamacterilstic time of the chemical reaction, 
E is the activation energy, R is the universal gas constant, 
T is absolute temperature and ¢ is the concentration of oxygen. 
The reaction rate was determined by utilizing an empirical 
expression developed by Frank-Kamenetski [Ref. 7]. 

The three field equations described above are presented 
below: 


(hy) Fiber Heat Transfer Equation 


5 A h,* aT. 
axl tk, i543 Cisne es Nes AHr(T,, b)-0 Case Cr) 
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(2) Internal Flow Heat Transfer Equation 








oT aT eZ dT 
3 a  - ao hes ; : a 
3x | Ka op oo i eeelicamon ee! 
eo Oxygen Transport Equation 
a (p29) ee. ok a6 (11.13] 





| T 
- — __D _ 4 , = L = 
ax BSx “Dp 3x aT, | 5x)? eC >) 3t 


Several of the parameters included in equations II.11, II.12, 
ome 11.135 are functions of temperature (i.e. thermal conduc- 
tivity, reaction rate). This temperature dependence was 
retained in the mathematical model and the parameters were 


changed with temperature. 


C. BOUNDARY CONDITIONS 

In order to make the mathematical model results as mean- 
ingful as possible, it was desired that boundary conditions 
De imposed that simulated the experiments of Fontenot [Ref. 1]. 


Vatikiotis imposed the following boundary conditions: 





3T 
(Kg + Kr) == “Gee Ee) ceqt -T *) at x = 0 (11.14) 
aul 
eee = - aT, CTL) C- «ose (7,* - TT") at x = Lb (11.15) 
ee st agi 0 (Ginleaies 
ai—_ OT 
 * 2 ae x = L (Ll 
o> —— - — 
Ba oo ome) cat x =0 Glee 6 | 
[a z 
- Qi "at x LE es ae og 
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The above boundary conditions provide a reasonable model for 
the one-dimensional problem. Boundary conditions II.14 and 
Mmieeies cONSider convection and radiation heat transfer from the 
Merous Graphite structure. Boundary condition II.16 sets the 
air entering the porous plate at ambient temperature. Boundary 
@ematttzons [1.18 and [1.19 are Dankwert conditions [Ref. 8] 


for flow through a porous media. 


D. ADDITIONAL ASSUMPTIONS 

The mathematical model is further based on two additional 
assumptions not previously described. 

(1) The temperature across each individual fiber is as- 
sumed to be constant. This assumption 1S quite reasonable due 
to “the small diameter of the gravhite fibers. 

(2) The temperature of the air in each individual pore 
is assumed to be constant. 

These two assumptions should not be interpreted as limiting 
fee temperature from fiber to fiber, or from pore to pore. 
The model allows temperatures to vary according to the equations 


described above. 


Ee weletHOD OF NUMERICAL SOLUTION 

The field equations form a system of three nonlinear, 
coupled partial differential equations. A Galerkin formulation 
of the finite element method is utilized to obtain solutions 
Ge Gie ftela equations. Linear interpolation fuctions are used 


ae ume basis functions. 
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The plate thickness (L}) is considered to form a closed 
domain (0, L) that is partitioned into (n-1) elements. This 
Brecess results in an n nodal point model. In accordance with 
the Galerkin method, a system of ordinary differential equa- 
eeenms 1s finally obtained. 

The resulting system contains 5n ordinary differential 
equations (graphite temperature, air temperature, and oxygen 
concentration at each nodal point). These equations retain 
[me Character of the original set of partial differential 
equations. 

Initial tests were conducted using a 21 nodal point model 
Meo-d.c.). Several instances of numerical instability were 
discovered for twenty one nodal point computer runs. Thirty- 
one nodal point and thirty-five nodal point models were tested, 
With convergence obtained within five percent. A thirty-one 
modal point model (93 o.d.e.) was selected for all subsequent 


mest runs. 
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Peele ian: FORMULATION 


It was desired that the thermal response of a composite 
meece De thoroughly analyzed. In order to carry out this 
objective, a test plan was formulated for use with the one- 
dimensional combustion model described in Section II. Initial 
efforts in formulating the test plan concerned selection of 


input parameters for variation during testing. 


A. EXTERIOR WIND VELOCITY (U,) 

)eceri1or air flow rate (U_) was immediately selected as 
a test parameter because of its obvious effects on average pore 
velocity and because it acts as a source of oxvgen for the 
combustion process. It was felt that relative wind velocities 
between ten and ninety knots would be representative of the 
velocities commonly experienced on aircraft carrier flight 
decks. Wind velocities of 10, 40, 70, and 90 knots were se- 


lected for testing. 


foe FRLeE THIGKNESS (L) 
Peatemente mess (lL) is a Critical parameter because it 
also effects pore velocity. <A consideration of Darcv's Law 


POunts out the effects of plate thickness on pore velocity, 


— Cini a 


Dp 


where fo is the pressure differential across the plate, m is 


<|s 
dé 


the specific permeability of the medium, and u is the dynamic 


oe 





[eeeeotey. NS Can be seen in equation III.1, pore velocity 

(.,) decreases hyperbolically as plate thickness increases. 

Mie efrect of plate thickness on pore velocity was felt to be 
Sete relevant because it could produce extremely high or ex- 
@eemely low pore velocities, both of which could result in 
@emeeeton. Plate thicknesses of 1.0 inch, 0.5 inch, 0.25 inch, 
and 0.125 inch were selected for testing. t was felt that 
these thicknesses were typical of composite plates under con- 


Sideration for use as airframe material. 


mee ANTERNAL PLATE GEOMETRY 

The pore velocity equation(I.1) indicated that both per- 
meability and porosity of the porous medium were significant 
factors in analyzing the thermal response of the composite 
plate. Porosity is a measure of void space and is indicative 
of the space available for oxygen, while permeabilitv represents 
miesability of a porous medium to allow flow. Equations [1.2 
and II.3 indicated that filament diameter ({d) and plv thickness 
(D) significantlv effect both porosity and permeability (see 
meere 1), ‘Filament diameters of .0002 inch, .0004 inch, .002 
inch, and .004 inch were selected for testing. Ply thicknesses 
meevodo5 mmen and .005 inch were utilized in the testing program. 
Mimese selections of filament diameter {d) and ply thickness (D) 
allowed for consideration not only of the effects of fiber 
diameter, porosity, and permeability, but also the effects of 


changing permeability while holding porosity constant. 


l@ wa 


suemese desimed to consider the effects of changing 


Merosity while holding permeability constant. In order to 


accomplis 
piy thick 


the same 


h this objective, specific filament diameters and 
messes were selected such that permeability remained 


in each case. The following selections accomplished 


ms task: 


ie d 
2 dad 
ee od 
The final 


matrix of 
ent cases 
-), air 
BP rawine 
time and 
(7) 
(2) 
(3) 
(4) 
(3) 
(6) 
oy) 
(8) 
(9) 
Therefore 
changes i 


System pa 


MOOS ctmen:,. D FUT eT Boverlal 


wOUGmeemen: DPD = 17007 inch 

——  U0Geeinen; D = .0007 inch 

S2eOwee 1S Dresemted as Table I]. This (4 x 28) 
computer runs resulted in production of 112 differ- 
for analysis. In addition to graphite temperature 
temperature (T.), and oxygen concentration, the 
System parameters were calculated as functions of 

position for each case: 

Filament diameter 

Porosity 

Permeability 

Pore velocity 

Internal Reynolds number 

Internal heat transfer coefficient 

Reaction rate 

Filament conductivity 

Deeruscton coefficient of oxygen 

, each of the computer runs produced data concerning 

n temperature and oxygen concentration and changes in 
rameters. 
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| Tapes Lil: TEST PLAN 


PLATE FILAMENT PLY CASE CASE CASE CASE 
THICKNESS DIAMETER THICKNESS A B C D 
fem [IN] [IN] [IN] [KNOTS] [KNOTS] [KNOTS] [KNOTS] 
1 0.5 004 005 U,=90 U=70 U,=40 U_=10 
2 1.0 004 005 U.=90 U_=70 U_=40 U_=10 
3 0.25 004 005 U.=90 U_=70 U_=40 U_=10 
4 0.125 004 005 U.=90 U_=70 U_=40 U_=10 
pS 0.5 002 005 U.=90 U.=70 U_=40 U_=10 
6 1.0 002 005 U.=90 U.=70 U_=40 U_=10 
7 0.25 002 005 U.=90 U.=70 U_=40 U_=10 
g 12s 002 005 U.=90 U.=70 U_=40 U_=10 
9 0.5 0002 0005 U.=90 U_=70 U_=40 U_=10 
10 ine 0002 0005 U.=90 U_=70 U.=40 U_=10 
_ 0.25 0002 0005 U.=90 U.=70 U_=40 U_=10 
Mm 860.125 0002 0005 U.=90 U.=70 U_=40 U_=10 
mm 0.5 0004 0005 U.=90 U.=70 U_=40 U_=10 
m §61.0 0004 0005 U.=90 U_=70 U_=40 U_=10 
m 0.25 0004 0005 U-=90 U.=70 U.=40 U_=10 
16 = 0.125 0004 0005 U-=90 U_=70 U_=40 U_=10 
m 0.5 003 004 U.=90 U_=70 U_=40 U_=10 
/ 1.0 005 004 U.=90 U_=70 U_=40 U_=10 
19 0.25 003 004 U.=90 U.=70 U_=40 U_=10 
Mm 86s 00.125 003 004 U.=90 U.=70 U_=40 U_=10 
1 ie 006 007 U.=90 U.=70 U_=40 U_=10 
2s. 0 006 007 U.=90 U_=70 U_=40 U_=10 
860. 25 006 007 U.=90 U_=70 U,=40 U,=10 
oS 0.125 006 007 ee 0 eA One) 0 
im o0.5 0002 .0007 U.=90 U_=70 U_=40 U_=10 
fh «1.0 0002 0007 U.=90 U_=70 U_=40 U_=10 
aw 0S 0002 0007 U.=90 U_=70 U,=40 U_=10 
om 0.125 0002 0007 U.=90 U_=70 U_=40 U_=10 
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ieee eeoUll so AND OBSERVATIONS 


In analyzing each computer run, particular attention was 
directed at the combustion-extinction phenomenon. Initial 
muemut imdneated that a partacular combustion flash point did 
Meteexist for prescribed initial conditions. Instead, it was 
surmised that a combustion continuum exists for a composite 
plate, ranging from immediate extinction, to immediate total 
mlamwe COmbUStion. This conjecture resulted in an additional 
research objective of attempting to determine the nature of 
this continuum. Placement of the different computer runs in 
the speculated combustion continuum was a necessarily subjec- 
tive process. The following general categories were utilized: 

(1) Immediate extinction 

oe) Slit cembustion, fast cooling 

06 6Seitcombustizon, slow cooling 

(4) Moderate combustion 

jy) Extensive combustion 

(6) Rapid combustion 

The effect of each of the input parameters on the combus- 
tion-extinction process was studied extensively. It was 
expected that by analvzing these effects for numerous computer 
runs, that the combustion continuum for a particular plate and 
prescribed initial conditions would be understood, and that the 


effects of svstem parameters would be discovered. This 
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information could then be utilized to select those design para- 
meters that improve the flame resistance of composite materials. 
Improved flame resistance would increase the survivability of 

a composite structure. 

A graphics package was inserted in the program to allow for 
ease of interpretation of results. Graphite temperature, air 
memperature, Oxygen concentration, reaction rate, and the 
Various system parameters listed above were all included in 
the plotting procedure. The plots are presented as functions 
Be plate location for particular times. 

The initial conditions prescribed for the 112 computer runs 
are presented in Table III. 

An initial graphite temperature of 1050°F was selected as 
a result of several preliminary test runs. These test runs 
indicated an ignition temperature in the vicinity of 1050°F 
for filament diameters (d) and plv thicknesses (D) of the same 
order of magnitude as those called for by the test plan. 

Initial air temperatures were selected to simulate a steady 
increase in temperature from ambient (50°F) at the entrance to 
MOSO-F in the plate interior. 

It was observed that initial oxygen concentration had 
little effect on the solution. The mathematical model imme- 
diately settled on the oxygen concentration appropriate for 
Pie Barticuldr test case, no matter which initial oxygen con- 


Ceomeration was utilized. 
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tect: INP EAL CONDITIONS 


GRAPHITE AGB 





NODAL OXYGEN 
POINT ee TEMPERATURE TEE ERATE | CONCENTRATION 
il Or 0 1050.0 505.0 .004 
¢ 0295.5 55 IGS 0 a 6 ho0e 0 .004 
5 0.06667 Oso. 0 25,050 .004 
4 0.10000 O50 0 6510500 .004 
5 Oe ie kegs Bo50 0 450.0 .004 
6 OR levers HOS. 0 55020 mone re! 
7 0.20000 Os 05, 0 650.0 .004 
8 O25) 5 5 ny S50 . 0 75020 .004 
9 2000 / POS: 0 850.0 .004 
10 U5 000.0 1050.0 950.0 .004 
1} 0] 34555 FOS 0'. 0 TO S070 .004 
ii? 0.36667 FOS0' . 0 OS 0 er0 .004 
iS 0.40000 HOS 0 OVS) 0 .004 
14 0.43555 O50 0 OS0r20 .004 
) a 0.46667 WOr5 0,. 0 J S06 .904 
16 0.50000 05 0) 20 HOS 070 .004 
a7 Op S555 1050-0 O50. 0 .004 
18 03,5666 / 5: 0.-0 hos 0e.0 .004 
19 0.60000 P50 20 Je S8) 0, .004 
20 0.65555 £@0) 0 JO) 3.0) 410) .004 
Zi 0.66667 10.50). 0 OS 0 0 .004 
22 0.70000 1510) SaO Lage 105 02..0 .004 
25 Oe? S55 1050.0 EOS.0.2-0 .004 
24 De 1066 / 05 .0..:0 Op 0s, 0 .004 
cas 0.80000 Oise U Fos.0". 0 .004 
26 55,555 £05,040 es. 0 .004 
a7 Oo6 0 6 / 0350-0 TORS 0; 0 .004 
28 0.90000 1050.0 GES Ol 10 .004 
29 CS eae BOs0 20 Oss O80 Ore 
30 0.96667 Ors 0-0 BOS07 0 Ove 
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Mieoec  ingwtial conditions produced results that ranged 
from immediate extinction to moderate combustion for the 112 
test runs. The influence of the input parameters on the com- 
bustion-extinction phenomenon was quite pronounced in several 


cases, these effects will now be considered. 


eee tS OF PLATE THICKNESS 

The test plan provided four plate thicknesses in order to 
explore the effects of this parameter on thermal response. 
The thicknesses were as follows: 

Sos Ul.125 inch 

ee 0-250 inch 

fed. 500 8 inch 

oe 1.000 inch 

It was generally observed that varying plate thickness pro- 
duced two extinction regions, existing on both sides of a 
combustion region. Slight to extensive combustion was observed 
Mor the test runs when in the combustion region. The two ex- 
tinction regions are the result of two different phenomena: 

meee Ghokine' of oxygen flow - pore velocity and oxygen 

Gemeentration are depleted to such an extent that 

combustion cannot be supported; 
(2) 'Blowout' - oxygen concentration and pore velocity 

are of such magnitude that convective heat transfer 


overcomes heat generation and cooling occurs. 
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Pmomewpeccaurie Crample of plate thickness effects consider 
test sets 1-4. The following parameters characterized the plate 
Geometry for these test sets: 

(1) Filament diameter = .004 inch 

Peewee thackness = .005 inch 
Tae. 25 inch plates (test set 4) experienced extinction for all 
meamemrOn Velocities. Figure 2 depicts the extinction process 
memecge OL125 inch plate, with an exterior wind velocity of 90 
Knots. Graphite temperatures are decreasing at all locations 
and oxygen concentration is high at all locations, ranging from 
.007-.008 lbm/FT*. These thin plates are characterized by high 
pore velocities, typically greater than 10,000 FT/HR. (see Figure 
5). These high pore velocities and oxygen concentrations indicate 
tmat extinction occurs as a result of the ‘blowout’ phenomenon. 
The high pore velocities enhance convective heat transfer pro- 
meme cme COOlimg effect depicted in Figure 2. The high pore 
velocities were expected for the thin plates because pore velo- 
Srty ws governed by Darcy's Law. As previously pointed out, pore 
velocity increases hyperbolically as plate thickness decreases. 

The .25 inch plates (test set 3) experienced moderate 
combustion for all velocities tested, with the exception of 
Beommentots. Figure 4, 5, and 6 depict a typical combustion 
eomemece 10r plates of this thickness. As these figures 
indicate, the combustion process initiates in a region removed 
from the plate entrance, eventually moving to the surface of 
eiemolete. Lt appears that some initial cooling of the graphite 


fibers is occurring in the plate entrance regions. The 
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eeeeuation is Only temporary though, until self-sustaining com- 
DieectOm Initiates in the plate interior. The oxygen supply 
follows the combustion sequence with all oxygen being consumed 
in the combustion region and zero oxygen concentration existing 
Mevyeond this region (see Figure 7). The oxygen concentration 
results explain the movement of the combustion region from 
the plate interior to the plate surface. The combustion pro- 
M@esoeiS Utilizing all the available oxygen in order to sustain 
itself, but additional oxygen could be used to enhance combus- 
tion. Therefore, the combustion region moves toward the 
Emerance, Where the most abundant oxygen supply exists. 

Peme VelOcities for the 0.250 inch plates are greatly 
memeca trom those of the previously discussed thinner plates. 
Velocities of 7000 FT/HR are typical. The increased thickness 
Peduced the pore velocity to such an extent that convecting 
heat transfer is no longer removing heat at a rate great enough 
to provide cooling, but the velocity is of sufficient magnitude 
—™@ Drovide enough oxygen to support moderate combustion in the 
Beer half of the plate. It is thought that combustion initiates 
in the plate interior and not on the plate surface because 
Gene velocity is high enough to provide some initial cooling 
at the plate entrance. It is further thought that continued 
reductions in pore velocity would result in reduction of the 
Somlee entrance region. A critical pore velocity could be 
reached where no entrance cooling is achieved, but the velocity 


me sWieticient to support combustion throughout the entire 
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Meee. POre velocity reductions beyond this critical point 
would result in insufficient oxygen supply to support combus- 
mon throughout the entire plate. Combustion initiation in 
Mme <tront half of the plate would characterize this situation. 
Meee results obtained from the .5 inch plates appear to support 
this speculation. 

The .5 inch plates (test set 1) experienced moderate com- 
bustion similar to the .25 inch plates for all velocities 
Mested, With the exception of ten knots. The combustion pro- 
Bess tOr these plates differed from that of the .25 inch plates 
in the location of the combustion initiation region. As can 
be seen in Figure 8, combustion initiates in the plate entrance 
region. Pore velocities of 3500 FT/HR were typical for these 
mmtees. itt appears that the reduction in pore velocity greatly 
meauecd the cooling effect in the plate entrance regions noted 
mummies. 2s inch plates, but that the velocity was still high 
enough to support combustion in this region. 

Miewleorinen plates (test set 2) experienced slight local 
combustion in the plate entrance region, with eventual 
eee tton Occurring for all velocities tested. In all cases, 
pore velocities were not sufficient to support combustion for 
an extended period of time. As can be seen in Figure 9, oxygen 
concentration is rapidly falling even with an exterior wind 
velocity of ninety knots. This lack of oxygen within the 
plate interior effectively chokes the combustion process 


BeGuitime if extinction at all velocities tested. 
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Mime errects Of thickness on thermal response are quite 
eyr@ent trom the above tests. Thin plates are desirable in 
that extinction by 'blowout' is frequently achieved, but thin 
Meaees are Still susceptible to combustion if average pore 
Mevocity is reduced to such an extent that the thermal response 
moves from the 'blowout' region to the combustion region. Re- 
@eece exterior wind velocity, or a change in internal plate 
femmetry that results in reduced permeability could both move 
mre tiermal response for a thin plate to the combustion region. 

Mieek plates are advantageous in that extremely low pore 
memocteites™are trequently achieved. This situation results in 
Peoeieremme “ot oxvgen flow that will lead to extinction. The 
eGhoking phenomena is not particularly desirable in that the 
pemec Interior does not cool as rapidly as when extinction 
Seems Dy “"blowowt'. High temperatures remain in the center 
Sememe plate for long periods of time. This type of situation 
probably occurred when 'smoldering' was observed in several of 
the tests conducted by Fontenot [Ref. 2]. 

Marek plates are also susceptible to combustion if high 
meeervor velocities occur or if plate geometry is altered so 
that permeability is increased. 

Ineanaly2ins the tests for various thicknesses, a good 
deal was learned about combustion initiation. It was learned 
Sree eve™inttiation région is controlled by pore velocity, 
Whteh is inflwenced by plate thickness. It was also learned 


emet tMe size of the combustion initiation region is quite 
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mem@citive to pore velocity, with a critical velocity existing 
Berea Particular plate geometry, such that the entire plate 


Mmeisexperience combustion if this velocity is present. 


fee EREECTS OF FILAMENT DIAMETER 

Altering fiber diameter, while holding all other system 
parameters Constant, again results in the three characteristic 
Mremmal response regions previously described (choking - com- 
bustion-blowout). Varying only fiber diameter results ina 
Smanege in the internal geometry of the plate. A decrease in 
fiber diameter (d) will effect an increase in both porosity 
and permeability, resulting in more space available for oxygen 
within the porous graphite structure and less internal resis- 
lemme tO Oxygen flow. This situation was explored in more 
detail in test sets 5-8. In these test sets fiber diameter 
was reduced from .004 inch (test sets 1-4) to .002 inch (test 
Sets 5-38). 

Eee tion by 'blowout' was observed for all thicknesses 
meemed at both 90 knots and 70 knots (Cases 5A, 5B, 6A, OB, 
74h, 7B, 8A, and 8B). Figure 10 (Case 5A) is an example of the 
thermal response observed for these test cases. The reduction 
in fiber diameter resulted in more oxygen per unit volume and 
Meomer pore velocities than those observed for similar cases 
Mmmeeest sets 1-4. This situation enhances convective heat 
Eranster producing the rapid cooling depicted in Figure 10 


(Case SA). Reduction of fiber diameter had a pronounced effect 
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mememe 0.25 inch plate and the 0.50 inch plate at both 90 
Knots and 70 knots. When fiber diameter was .004 inch, both 
plates experienced moderate combustion. When fiber diameter 
was reduced to .002 inch the thermal response was altered from 
the combustion region to the extinction by 'blowout' region. 

Moderate to extensive combustion was observed in the 
following test cases: 

wees 26 (0 =40, L=0.50, d=.002, D=.005) 

emmeaorcom oe (U=10, L=0.50, d=.002, D=.005) 

em eoP (UU =10, L=1.0, d=.002, D=.005) 

eemeeetere 7D (Ue=10, L=0.25, d=.002, D=.005) 
Pugeme lil (Case 5C) depicts an extensive combustion situation, 
Mauic Figure 12 (Case 5D) presents a more moderate thermal re- 
Sponse. In these situations it appears that the higher porosity 
and permeability is enhancing the combustion process by pro- 
Viding more oxygen. The external velocities and plate thick- 
messes present in these cases result in pore velocities that 
Support combustion, but are not of sufficient magnitude so as 
to provide a cooling effect. 

Extinction was again observed in the remaining test cases: 

Oe peteee OC (U =40, L=0.25, d=.002, D=.005) 

Poeeeeee soe (U0 -40, L=0.125, d=.002, D=.005) 

Pepeeee SP (U_=-10, L=0.125, d=.002, D=.005) 
In these cases very high pore velocities (3500 FT/HR -55,000 
me jana Oxysen concentrations were present. The extinction 


mechanism was once again ‘blowout’. It should be noted that 
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m@esecond Characteristic extinction mechanism, 'choking', due 
to insufficient oxygen to sustain combustion, was not observed 
tor these test cases. 

Baummes Speculated that the ‘choking' extinction region 
could be observed by conducting additional tests at greater 
thicknesses, thereby decreasing the pressure differential 
across the plate. The reduced pressure differential would pro- 
duce lower average pore velocities. In order to confirm this 
speculation additional tests were conducted at the following 
@nicknesses: 

(1) 2 inches 

ie > inches 

ji) %S inches 
Filament diamter and ply thickness remained at .002 inch and 
.005 inch respectively for these tests and external velocity 
was 90 knots. 

Mee 2 inch plate resulted in extensive combustion. Extinc- 
tion by ‘blowout’ had been previously observed for all thick- 
Messes tested at 90 knots. This result for the two inch plate 
confirmed the expected combustion region. Combustion initiation 
Ee@emtic two inch plate occurred in the center region of the 
plate, with a cooling effect noted in the vicinity of the plate 
entrance. The location of the combustion initiation region 
indicates that the pore velocity is slightly higher than the 
velocity necessary to support combustion throughout the entire 


pate . 





The five inch plate resulted in extensive combustion once 
again, but combustion initiation now occurred near the plate 
Surtace. This change in position of the combustion initiation 
region indicates that pore velocity has been reduced below 
Mme critical point where combustion would occur throughout the 
entire plate. These results correlate well with the results 
meeeviously discussed in section IV A. 

Mm@enwo inch plate resulted in an extinction situation as a 
result of a ‘choking' of oxygen flow. As depicted in Figure 
13, oxygen concentration rapidly decreases. Slight combustion 
occurs near the plate entrance, but the oxygen flow is not 
sufficient to support even this slight combustion for an ex- 
tended period of time. 

Tire three additional thickness tests described above con- 
firm the existence of two extinction regions and one ignition 
region for the more porous plate, as previously speculated. 

Tests sets 5-8 indicate that decreasing filament diameter 
mnile holding other parameters constant results in an increase 
mie both permeability and porosity. This situation typically 
results in either extensive combustion or extinction by 'blow- 
ome’. as described above, for plates with thicknesses of less 
Maan 1 inch. These results would indicate that permeability 
gad porosity significantly effect the internal flow 
characteristics thereby controlling the thermal response of 
the composite plate. As previously noted, pore velocity is a 


function of both permeability (m) and porosity (p). 
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Meer rPeClS OF PERMEABILITY 
omoseer tO analyze the effects of permeability, test sets 
meebo are analyzed in comparison with test sets 1-4. The fol- 


lowing information contrasts these two test groups: 


eeeeceety. COMPARISON OF TESTS SETS 1-4 WITH TEST SETS 13-16 


PARAMETERS TEST SETS 1-4 TEST SETS 13-16 
Filament Diameter mO Oae pre 

Ply Thickness 00s" .0005" 
Permeability 20 Bent) Meise sis) se NiO he 
Porosity ~497 -497 





As can be seen from the above information, porosity remains 
constant while permeability decreases by a factor of 100 for 
test sets 13-16. It was expected that this significant reduc- 
peorein permeability would aiter the thermal response of the 
composite plate. 

Me "eemoust ion was “observed for any of the tests conducted 
mretest sets 135-16. Extremely low pore velocities were charac- 
Mermstic Of all the test cases. Pore velocities ranged from 
Seer iyHR for the 1.0 inch plate with an exterior velocity of 
Memicots, to i134 FT/HR for a .125 inch plate with an exterior 
Mind velocity of 90 knots. These velocities are not of 
sufficient magnitude to support even slight combustion. 


Extinction occurs in all cases as a result of ‘choking' of the 
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Meee low. The 1 inch, 0.5 inch, and 0.25 inch plates, all 
Smperienced zero oxygen concentration within the plate interior. 
Mere complete depletion of oxygen occurred at the locations 


indicated in Table V. 


TABLE V; OKYGEN DEPLETION LOCATIONS 


Pes THICKNESS K/L-ZERO OXYGEN 





Table V indicates that for plates greater than 0.25 inch 
not only is pore velocity insufficient to support combustion, 
Bee permeability and pressure differential across the plate 
are such that oxygen does not even flow across the entire 
plate. This is a very desirable situation from the standpoint 
Gr Combustion prevention. 

A comparison of test sets 1-4 (d=.004 inch, D=.005 inch) 
With test sets 13-16 (d=.0004 inch, D=.0005 inch) shows the 
significance of permeability to the combustion process. In 
mest sets 1-4, combustion was observed in nine of the sixteen 
tests conducted. In test sets 13-16, combustion was not ob- 
merved in any of the tests conducted, despite the fact that 
porosity remained the same. These tests indicate that a mini- 


mum porosity is a necessary condition for plate combustion, 








Meme It 1S not a sufficient condition. Porosity is an indica- 
tor of space available for oxygen, while permeability is a 
Measure of the ability of a porous medium to allow flow (i.e. 
high permeability implies low resistance to flow). Previous 
meees Mave indicated the critical nature of pore velocity in 
its relationship to the combustion process. Therefore, per- 
meability will have a significant impact on the resulting 
thermal response of the composite plate. 

Mestescts 15-16 (d=.0004 inch, D=.0005 inch) also provided 
some insight into the cooling process of the composite plate. 
Figures 14, 15, 16, and 17 depict graphite temperatures for 
meee anen, O.5 inch, 0.25 inch, and 0.125 inch plates re- 
mmeectively. these figures indicate that plate cooling procedes 
more quickly for the thinner plates. A consideration of the 
pore velocities and the oxygen concentrations that characterize 
Meese test cases is relevant to the cooling process. Table VI 
mmesents this information for the four plate thicknesses tested. 
The oxygen concentrations presented in Table VI all occurred 
Mery tvapidly after the problem start time (t less than 1.0 
second). 

Teele VI illustrates pore velocity increases as plate 
mamckmess decreases. This situation occurs because as plate 
thickness decreases the pressure gradient across the plate 
increases. The increased pressure gradient results in increased 
Dame velocity. Despite the fact that pore velocity is increasing, 


it is not sufficient to support combustion, even for the 0.125 
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Meermomace. Pore velocities of at least 500 FT/HR character- 
ized previous test cases that experienced combustion. As can 
be seen from Table VI, the pore velocities for the cases under 


consideration are well below this value. 


ae oe Vl 


mere VELOCITY AND OXYGEN CONCENTRATION (TEST SETS 13-16) 


PLATE SIZE PORE@ yr cO.G.l 1 ¥ Oe GEN CONGENERAT LON 
mao. inch i ISIE eI OR ata) lhe" 0 Sos 
feo inch So Fly ae ede = B0oo5 > 
0.25 anech O6..0 VEHeAnk ) Sate /L = 20C07 
eto inch Bot. 0 Rr HR .005 lbm/FT? 
(throughout entire 
plate) 


Even though the pore velocities experienced for these test 
cases were of insufficient magnitude to support combustion, it 
appears that they were high enough to provide some cooling to 
Mie plate through convective heat transfer. In fact, the 
Booling was quite significant for the 0.25 inch plate and the 
0.125 inch plate as Figures 17 and 18 demonstrate. 

Bematpoe@rs that ‘choking’ of oxygen flow and rapid cooling 
Secur Simultaneously for thin plates with reduced permeability. 
The reduced permeability provides the mechanism for preventing 
Mapeial combustion by reducing pore velocity to such a level 
that combustion cannot be sustained. The thin plate thickness 


results in a pore velocity that is sufficiently high to cool 
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mime Gbhate through convective heat transfer. These results in- 
meeaee that the composite plate designer can not only design 
mr combustion prevention, but also for desired cooling rate. 
This can be accomplished by not only designing for a pore 
velocity such that the thermal response is in the extinction 
by choking region, but also designing the plate geometry such 
that the oxygen flow is not completely choked within the plate 
MmectriOr. it is believed that a threshold pore velocity exists 
for any particular composite plate, such that combustion will 
not be supported and a maximum cooling rate will be achieved. 

Mmmeecotoesets 9-12 (d=.0002, D=.0005) the vermeability 
Mech@entliy hieher (.625 x al er) than that for test sets 
13-16 (.29 x iit aig ET). This increase in permeability 
Mmecultod@ in Nigher pore velocities in all test cases. This 
Situation resulted in moving the thermal response from the 
extinction by ‘choking’ region to the combustion region for 
meveral test cases. Combustion was observed in eleven of the 
mexeemm tests Conducted for these test sets. 

foal the test cases in which combustion occurred, ini- 
Deeeenmetaok place in the vicinity of the plate entrance. It 
was also observed that the combustion process was not of a 
moDUSit Mature, with slight to moderate combustion observed in 
Srenwedse. Fisure 18 (case 9C) presents the thermal response 
meereosoeemcha plate, with an exterior velocity of 40 knots, 
Seemcmcypical of the combustion process observed for these 


test sets. Additionally, oxygen concentration was rapidly 
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Mepetead in the vicinity of the plate entrance for all con- 
PiisceroOn Situations. Figure 19 (Case 9C) depicts the oxygen 
Mmepletion process observed in combustion situations. The 
combustion initiation region and the rapid oxygen depletion 
both indicate that the combustion process is immediately 
utilizing all oxygen available. These results indicate that 
the permeability is retarding the thermal response. A further 
increase in permeability would result in less resistance to 
oxygen flow, which would produce pore velocities of greater 
magnitude. The higher pore velocities would move the thermal 
response further into the combustion region, resulting in a 
more extensive combustion situation. 

It should be noted that porosity was substantially higher 
Mrecest sets 9-12 (p=.9745) than in test sets 1-4 (p=.497) and 
test sets 135-16 (p=.497). Previous tests indicated that the 
thermal response was not as sensitive to changes in porosity 
as it was to changes in permeability. In these tests though, 
permeability was altered with porosity. In order to completely 
mies tiadte the effects of porosity, test sets 1/-28 were con- 
ducted. These tests were designed to vary porosity, while 


Pomemne permeability and all other input parameters constant. 


we) ERPRECTS OF POROSITY 

Mestesets 17-28 investigated the effects of varying poro- 
euley, While holding all other parameters constant. Table VII 
Summarizes the parameters that characterize the various test 


Bets. 











imeeeee vil: 6©6PRREMETERS FOR TEST SETS 17-28 


FILAMENT PLY a 
TEST SETS i eT PERMEABILITY | POROSIT 
| 9 z 


008 inch S00 4aanch 
S00 iGeeamch 00g ch 
rOCge2s ineh TELONC HE Saigielel 





Pretest oets 21-24 (p=.420) 

As Table VII indicates, porosity and permeability for 
test sets 21-24 are the same order of magnitude as those experi- 
enced in test sets 1-4 (m = .29 x gna. p = .497). As ex- 
pected, test sets 21-24 presented thermal responses very similar 
to those observed for test sets 1-4. Combustion was observed 
in nine of the sixteen tests conducted for these sets, with 
combustion initiation occurring in a region removed from the 
plate entrance. It was also observed that the combustion pro- 
cess was not of a robust nature, with slight to moderate 
combustion observed in each ignition case. Extinction as a 
result of the 'choking' of oxygen flow was observed for the 
eometmea 0.5 inch, and 0.25 inch plates when the exterior 
Weloeity Was ten knots. Extinction as a result of 'blowout' 
See ac4am observed for the 0.125 inch plates at all velocities 


tested. 
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eee oets 17-20 (p=.560) 

Meese sets 1/-20. porosityawas increased to .560, 
mmile permeability remained at .50 x io Fee This increase 
in porosity (approximately 350 percent) results in more void 
meereceavawlable for oxygen. Combustion occurred in cases 
Similar to those that experienced combustion in test sets 
meas. EXtinction was again observed for the 1.0 inch, 0.5 
miem, and 0.25 inch plates with exterior velocity equal to 
Mem Knots, and for the 0.125 inch plates at all velocities. 
The thermal response observed in these test sets, when com- 
bustion took place, was markedly different than that observed 
for the combustion cases in test sets 21-24. Extensive com- 
bustion was observed in several instances, especially for the 
0.25 inch and the 0.50 inch plates. Figures 20 (Case 174A) 
and 21 (Case 19A) are typical of the thermal response observed. 

It appears that the increase in porosity from .420 (test 
mets 21-24) to .560 (test sets 17-20) did not induce combustion 
mom@celr in any additional test cases, but it did impact on 
the thermal response. The combustion process is of a more 
mopmet nature. The additional space available for oxygen 
Mee orts the combustion process once it initiates. 

bwetescteoets 25-28 (p=.935) 

Mmeeestts sets 25-28, porosity 15 increased to .9535, 
indicating a great deal of void space available for oxygen. 
Somemstaen Was again observed in several instances. The sig- 


Migicant increase in porosity produced a dramatic change in 
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the speed of the observed thermal response. Figure 22 (Case 
25B) is indicative of the rapid thermal response observed 
when combustion took place in these high porosity tests, as 
mmeure a2 depicts, graphite temperatures in excess of 1250°F 
Bee achieved in less than 0.25 seconds. These tests again 
indicate that an increase in porosity does not induce combus- 
tion to occur in any additional test situations, but it does 
impact on the nature of the thermal response. High porosity 
(p greater than 0.90) appears to result in rapid, extensive 
Semoustion Situations if permeability is such that pore velo- 
city will support the combustion process and allow combustion 
to initiate. 

In summarv, it is thought that permeability, and its 
mesooctatead impact On pore velocity, determines whether the 
thermal response will be located in the combustion region or 
in one of the two extinction regions that exist on either side 
Of the combustion region (extinction by 'choking' of oxygen 
flow or extinction by 'blowout'). If the thermal response is 
Mecated im the combustion region, it is thought that poro- 
sity controls the speed and robustness of the combustion 


process. 


Peeeeeeeeeis OF EXTERIOR VELOCITY 
EPapenior velocity is the driving force for the internal 
oxygen flow that is required for both the combustion and the 


extinction processes. Even though exterior velocity provides 
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this driving force, neither extremely high nor extremely low 
MenOcitiles are sufficient to guarantee extinction. Previous 
meres stON nas indicated that pore velocity is a critical fac- 
tor in analyzing the thermal response. Permeability, plate 
BorekKmess, and porosity are all factors which can control and 
alter interior pore velocity to such an extent that the thermal 
response may be contrary to preconceived expectations based 
solely on the knowledge of exterior velocity. 

Exterior velocity does appear to effect the combustion 
initiation region. In several test cases, with high exterior 
velocities (1.e. 90 knots) initial combustion was observed in 
the latter half of the plate, while inthose cases with low 
exterior velocities, initial combustion occurred in the proxi- 
mity of the plate surface. These results indicate that even 
though plate geometry results in a pore velocity that supports 
combustion, the exterior velocity is controlling where the 
initial combustion will occur. High exterior velocities result 
mameeeeeenetlow that provides cooling at the plate entrance. 
Initial combustion then occurs in the latter plate regions after 
the oxygen temperature has reached a sufficient level so that 
heat generation overcomes convective heat transfer. Figure 5 
feweers COmpuUStion initiation in thelatter plate regions. 

Pewmextcerior velocites are not sufficient to produce the 
Mmmm coolans eftcts described above. If the pore velocities 
Presstiitieient to support combustion, ignition will commence 


in the proximity of the plate surface. 
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mmeewese SswtWations 1n Which extinction occurs as a result 
of 'blowout', higher cooling rates were observed with higher 


exterior velocites. This result was as expected. 


48 








V. CONCLUDING REMARKS 


Mimmecst results indicate that permeability, porosity, 
melament diameter, plate thickness, and exterior wind velocity 
fmeiave Significant effects on the combustion process. In 
attempting to understand the thermal response of a composite 
meoree, it must be realized that two extinction regions exist 
on either side of a combustion region. Figure 23, depicts 
this situation. For a given temperature, say T = 1050°F, ex- 
tinction will result for pore velocities less than Be Oper Ie@ie 
memew Velocities greater than uc. For u_A < u_ < uc, com- 

Pp P Pp Pp 
Bustion will result for all T > 1050°F. Thus for a given 
Memmperattre, and fixed design, there 1s extinction at low pore 
velocites (up < Basle followed by combustion for oe < Bs < wee 
gama again followed by extinction for a Bae: 

Pommmonre velocities less than Leen the system is air poor, 
and extinction occurs because heat generation is less than the 
heat transferred out of the system as a result of convection. 
As bs 1s increased beyond Oe heat generation exceeds heat 
transferred and combustion with increasing system temperature 
met Occur. The minimum ignition temperature occurs at eee 
Increasing a above oe Pesubeswmine Imercasing lonitlon»*tempera- 


meres. FOr pore velocities greater than uo extinction will 


Meer once again since heat transferred exceeds heat generated. 
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Memer@ompOrnt, the system is air rich. That is, there is an 
excess of oxygen above what is required to sustain combustion 
and ‘blowout’ will occur. 

Memes been established that an ignition curve similar to 
meet discussed above for Figure 25 exists for any composite 
meate. Average pore velocity is selected as the abscissa for 
Bme Curve because of the profound effects observed as a result 
of variation of this parameter during the test runs. Average 
Bore velocity is influenced by exterior wind velocity, plate 
thickness, and interior plate geometry. 

If a composite plate is designed to retard combustion by 
Semecime’ Of Oxygen flow, then high exterior velocities are 
not desirable because they may be of such magnitude that the 
thermal response is moved from the ‘'choking' extinction region 
to the combustion region. Similarly, if the designer intends 
combustion prevention to be achieved as a result of 'blowout', 
high exterior velocities would be desirable. [In short, the 
Meals of the composite material designer must be understood 
me cOmoustion is to be prevented and a desirable cooling rate 
mememe porows graphite structure is to be achieved. 

Additional, extensive testing is required in order to 
further quantify the thermal response of composite plates of 
Various geometries. Tests conducted with additional exterior 
Mand velocites and under different initial conditions are also 
memitered. &n exploration of the ignition process may prove 


enlightening. This task could be accomplished by beginning 
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Mith an extinction situation and gradually increasing initial 
Mmegatee temperature for each successive computer test run. 
Combustion initiation could be observed along with the robust- 


ness of the combustion process. 
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